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Medium modification of jet fragmentation function due to multiple scattering and in- 
duced gluon radiation leads directly to jet quenching or suppression of leading particle 
distribution from jet fragmentation. One can extract an effective total parton energy 
loss which can be related to the total transverse momentum broadening. For an expand- 
ing medium, both are shown to be sensitive to the initial parton density and formation 
time. Therefore, one can extract the initial parton density and formation time from si- 
multaneous measurements of parton energy loss and transverse momentum broadening. 
Implication of the recent experimental data on effects of detailed balance in parton energy 
loss is also discussed. 

1. Introduction 

In high-energy heavy-ion collisions, one of the pressing issues is the direct determination 
of the initial parton density of the dense medium and its formation time. Jet quench- 
ing or attenuation of leading particles from jet fragmentation can provide an effective 
tool to measure directly the parton density of the medium with which the jet interacts 
strongly during its propagation. Recent theoretical studies |H2|3|4|5] have shown that 
parton energy loss via induced radiation is directly related to the gluon density of the 
medium. The attenuation will suppress the final leading hadron distribution giving rise 
to modified parton fragmentation functions [6 . By measuring the medium modification 
of the fragmentation function one can thus extract the effective parton energy loss. Most 
importantly, one can compare the effective parton energy loss extracted from heavy-ion 
collisions to that of cold nuclei and extract the initial parton density of the hot dense 
medium relative to a cold nucleus. 



2. Parton Energy Loss in Cold Nuclei 

In deeply inelastic scattering (DIS) off a nuclear target, a quark suffers multiple scatter- 
ing with other nucleons inside the nucleus after it was knocked out of its parent nucleon. 
The induced gluon bremsstrahlung before hadronization leads to suppression of leading 
hadrons and gives rise to a modified quark fragmentation function. Including the leading 
twist-4 contributions from double scattering processes, the nuclear medium induced cor- 
rection to the fragmentation function is found to be similar to the radiative correction in 
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the vacuum except that the normal sphtting functions are replaced by the modified ones. 
These modified splitting functions are found to depend on the quark-gluon correlation 
function T^(a:, Here 

T^(x, xl) = dy^dy2e'^''^''^^P^y~+'''^P^^y^~y^\l - e'^^^^^^^ )(1 - g-^^^P^^^^-i'i")) 

\ {A\ij,{Q)^^ F^^{y,)F^'^{y^)Uy')\m-y2ny- -Vi) (1) 

and qA{x) is the twist-two quark distribution of the nucleus. Such a two-parton correlation 
function contains essentially four independent twist-4 parton matrix elements in a nucleus 
[xj. = {kl)/2p+q-z{l-z)]. The dipole-like form-factor (1 -e-^^^P^2^2")(l - e-*^^P^(^"-2^r)) 
arises from the interference between the final state radiation of the 7*g scattering and the 
gluon bremsstrahlung induced by the secondary quark-gluon scattering. By generalizing 
the factorization assumption to these twist-four parton matrices, we have 

T^g{x,XL)/qA{x) ^ C{Q^)mNRA{l - e'^i/^^), (2) 

with a Gaussian nuclear distribution p(r) ~ exp(— 
xa = ^/itinRa, and is the nucleon mass. 

Since the two interference terms in the dipole-like form-factor involve transferring mo- 
mentum xlp'^ between different nucleons inside a nucleus, they should be suppressed for 
large nuclear size or large momentum fraction xl. Notice that Tf = l/x^p^ is the gluon's 
formation time. Thus, xl/xa = La/tj, with La = Ra^n/p~^ being the nuclear size in the 
infinite momentum frame. The effective parton correlation and the induced gluon emis- 
sion vanishes when the formation time is much larger than the nuclear size, xl/xa <^ 1, 
because of the Landau-Pomeranchuck-Migdal (LPM) interference effect. Therefore, the 
LPM interference restricts the radiated gluon to have a minimum transverse momentum 
~ Q"^ I^nRa ~ Q^/A^/^. The nuclear corrections to the fragmentation function due 
to double parton scattering will then be in the order of a^A^/^/^y ~ ttgA^/^/Q^, which 
depends quadratically on the nuclear size. For large values of A and Q^, these corrections 
are leading; yet the requirement £^ <C for the logarithmic approximation in deriving 
the modified fragmentation function is still valid. 

Shown in Fig. ^ are the calculated nuclear modification factor of the fragmentation 
functions for ^'^A^ and ^'^Kr targets as compared to the recent HERMES data jH]. The 
predicted shape of the z- and and energy dependence agrees well with the experimental 
data. A remarkable feature of the prediction is the quadratic A^^^ nuclear size dependence, 
which is verified for the first time by an experiment. The only parameter in our calculation 
is found to be ^(Q^) = 0.0060 GeV^ with as(Q^) = 0.33 at ~ 3 GeV^. 

If one defines theoretically the quark energy loss as that carried by the radiated gluons, 
then the averaged total energy loss is, 

= u{^z,) ^ Ca'M)mMR\{CAlN,)?,\n{l/2xB). (3) 

With the determined value of C, {xb) ~ 0.124 in the HERMES experiment |8jj and the 
average distance {La) = Ra\J'^I't^ for the assumed Gaussian nuclear distribution, one gets 
the quark energy loss dE/dL f» 0.5 GeV/fm inside a Au nucleus. 
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Figure 1. Predicted nuclear modification of 
jet fragmentation function is compared to the 
HERMES data [Hj on ratios of hadron distri- 
butions between A and D targets in DIS. 



Figure 2. Calculated nuclear modification 
factor compared to PHENIX data [TUj . 
The lower dashed line used the energy loss 
as given in Eq. 0. 



One can also calculate the transverse momentum broadening of the quark jet which is 
also related to a similar twist-four parton matrix elements. 

= ^ C^mj^R^ ^ O.OllAV^GeV^ (4) 

Nc qA{xB) 

3. Initial Parton Density in Au + Au at RHIC 

To extend our study to jets in heavy-ion collisions, we assume ~ /x^ (the Debye 
screening mass) and a gluon density profile piy) = {tq/t)0{Ra — y)po for a 1-dimensional 
expanding system. Since the initial jet production rate is independent of the final gluon 
density which can be related to the parton-gluon scattering cross section SJ [ asXTG^xx) ~ 
fi'^CTg], one has then 



asTfl{xB, XL 



p dyagp{y)[l-cos{y/Tf)], (5) 



W"^^X„ rfrp(r)(r-ro)ln— (6) 



where Tf = 2Ez{l — z)/l\ is the gluon formation time. One can recover the form of 
energy loss in a thin plasma obtained in the opacity expansion approach 

_ -nCaCAal [Ra , ^, 2E 

TjJ? 

where we have assumed (Jg ~ Ca27ra^//i^ (Ca=l for qg and 9/4 for gg scattering). 

Neglecting the logarithmic dependence on r, the averaged energy loss in a 1-dimensional 
expanding system can be expressed as [9J 

dE dEo2To dEo TTCaCAal „ , 2E 
dL dL Ra dL 2 r^p^ 
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where dE^/dL is the energy loss in a static medium with the same gluon density po as in a 
1-d expanding system at an initial time tq. Because of the expansion, the averaged energy 
loss {dE/dL)id is suppressed as compared to the static case and does not depend linearly 
on the system size. Similarly, one have the total transverse momentum broadening in an 
expanding system, 

(Agi) = 27iCaal ^ drpir) = (Agl)o^ In ^; (Agi)o = 2nCaalpo (8) 
J TO JtLa To 

Fitting the PHENIX data [TU] on suppression of large ttq spectra in Fig. |21 yields 
{dE/dL)ict 0.341ni?/ln5 GeV/fm. Taking into account the expansion, this would be 
equivalent to {dE/dL)o = 0.3A{Ra/2to) In i?/ In 5 in a static system with the same gluon 
density as the initial value of the expanding system at tq. With Ra ~ 6 fm and assuming 
To ~ 0.2 fm, this would give {dE/dL)o ^ 7.3 GeV/fm for a 10-GeV parton. Since the 
parton energy loss is directly proportional to gluon density of the medium, this implies 
that the gluon density in the initial stage of Au + Au collisions at tq = 0.2 fm/c is about 15 
times higher than that inside a cold nucleus. To extract the parton density and the initial 
formation time, one has to measure the energy loss and qx broadening simultaneously. 

In Fig. 121 we also show the predicted suppression factor at large px with the same 
logarithmic energy- dependence of the energy loss as in the cold nuclear matter which 
gives a suppression factor that increases with px- However, new RHIC data [TU] at 200 
GeV show a constant suppression factor at larger pt- This indicates the importance of 
the detailed balance ^T] which gives much stronger energy dependence of the energy loss. 
We show as the lower dashed line in Fig. El the result for an effective parton energy loss 

AE (x{E/fi-1.6y-^y{7.5 + E/fi) (9) 

which is parameterized according to the result from Ref. JT] where both stimulated 
gluon emission and thermal absorption are included in the calculation of the total energy 
loss. The detailed balance between emission and absorption reduces the effective parton 
energy loss and on the other hand increases the energy dependence. The threshold is the 
consequence of gluon absorption that competes with radiation that effectively shuts off 
the energy loss. The parameter fi is set to be 1 GeV in the calculation shown in Fig. |21 
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